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^ ' Abstract 

m : 

Charm photoproduction rates off nuclei below the nucleon threshold are estimated using 
the phenomenologically known structure functions both for x < 1 and x > 1. The rates 
t-h ■ rapidly fall below the threshold from values ~10 pb for Pb close to the threshold (at 7.5 

GeV) to ~lpb at 6 GeV. 

m ■ 

■ 1 Introduction 

o : 

In view of the envisaged upgrade of the CEBAF facility up to 12 GeV it becomes important 
to have relatively secure predictions about the production rates of charm on nuclear targets 
below the threshold for the nucleon target. This note aims at such predictions. From the 
start it has to be recalled that the dynamical picture of charm production at energies close 
P-i. to the threshold is much more complicated than at high energies. On the one hand, the 

production rates cannot be described by the standard collinear factorization expression but 
involve gluon distributions both in x and transverse momentum, which are unknown at small 
transverse momenta (in the confinement region). On the other hand, in the immediate 
vicinity of the threshold the simple photon-gluon fusion mechanism of charm production 
becomes overshadowed by multiple gluon exchanges [1] and formation of colourless bound 
states with lower mass, as compared to open charm [2]. As a result, theoretical studies of the 
charm production near the threshold are known to be notoriously difficult. To avoid all these 
complications we study charm production at not too small distance from the thresholds, where 
all effects due to multiple gluon exchanges are hopefully small. To be more concrete, taking 
the charmed quark mass m c = 1.5 GeV, we have the threshold of open charm photoproduction 
on the proton target at the incident photon energy E{ h = 7.78 GeV and on the deuteron 
target at 5.39 GeV. As we shall argue (see Appendix), our treatment is hopefully valid at 
incident photon energies E excluding regions at distances ~ 0.3 GeV from these thresholds, 
say, 5.7 < E < 7.5 GeV. Closer to the thresholds multigluon exchanges and bound state 
formation may change our predictions considerably. 

In fact all estimates show that the ratio of the double to single gluon exchange contribu- 
tions is determined by the parameter a s (M 2 )[/i/(m c Ax)] 2 , where fi is the light quark mass, 
M is the mass scale for the gluon coupling to the charmed system and Ax is the distance 
from the threshold for the scaling variable of the produced charmed system ( [1], see also 
Appendix). Obviously as Ax — > multiple gluon exchanges begin to play the dominant role. 
They also bind the produced charmed quarks into colourless bound states, thus changing the 
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threshold value. However, with Ax growing their contribution becomes suppressed. Taking 
the effective fi ~ Aqcd = 0.3 GeV, and M = 2m c so that a s (M 2 ) ~ 0.3 we find the suppres- 
sion factor of the order 75(Ax) 2 . This factor determines the region where one can neglect 
all multiple gluon exchanges and bound state formation. In a more detail this suppression 
factor will be discussed in Appendix. 

Neglecting multiple gluon exchanges, subthreshold production will be determined by the 
gluon distribution in the "cumulative" region (x > 1 at large energies or massless target). As 
mentioned, in the low energy region the charm production rate involves the gluon distribution 
not only in x (collinear factorization) but also in both x and k\ (k± factorization). We shall 
use the simplest approximation about the structure of this combined distribution, assuming 
the dependence on x and k± factorized. For the latter dependence we shall use a monopole 
form, prompted by the perturbation theory. As for the former one, we shall exploit the 
existing (scarce) data on the nuclear structure functions at x > 1, from which we shall 
extract the relevant gluon distribution, using the DGLAP evolution equation. With all these 
simplifications, we hope to be able to predict the rates up to factor 2-=-3 1 . 

2 Kinematics and cross-sections 

2.1 The charm production cross-section 

Consider the exclusive process 

-f + A^CC + A*, (1) 

where A is the target nucleus of mass itia and A* is the recoil nuclear system of mass m* A . We 
denote the total mass of the CC system as M. Obviously M > 2m c where m c is the mass of 
the C-quark, which we take as 1.5 GeV. The inclusive cross-section for charm photoproduction 
is obtained after summing over all states of the recoil nuclear system. 

We choose a reference system in which the target nucleus with momentum Ap is at 
rest and the incoming photon with momentum q is moving along the z-axis in the opposite 
direction, so that q+ = q±_ = 0. The photoproduction cross-section corresponding to (1) is 
then obtained via the imaginary part of the diagram in Fig. 1 as 

a A ^ =Aj ^ 5({A P - k? - m A 2 )x ) ' o g {M\ k 2 ). (2) 

Here TAA*{k 2 ) is the vertex for gluon emission from the target; a g (M 2 ,k 2 ) is the photopro- 
duction cross-section off the virtual gluon of momentum k. We have also introduced the 
scaling variable for the gluon as x = k + /p + . Due to q + = this is also the scaling vari- 
able for the observed charm. Note that this definition, which is standard at large energies 
and produced masses, is not at all standard at moderate scales. In particular this x does 
not go to unity at the threshold for the nucleon target. Rather the limits for its variation 
converge to a common value 0.76. For the nuclear targets with A » 1 its minimal value 
at the nucleon threshold is well below unity (~ 0.64). One should have this in mind when 
associating this x with the gluon distribution: it follows that for a nuclear target, for energies 
going noticeably below the nucleon threshold, the cumulative region (prohibited for the free 
nucleons kinematics) includes not only values of x above unity but also a part of the region 
x < 1. 



1 Some crude estimates were earlier reported in [3]. Their order agrees with the present calculations at 
energies not too close to the threshold. 
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- A I m **$L* ^ ] ^M-X). (3) 



We use the 5-function to integrate over k- to obtain the cross-section (2) as 

-I } 2 

2{A-x){2itf ^ k 2 
In these variables we find 

M 2 = xsi + k 2 , si = 2pq (4) 

and 

k 2 = X Am 2 - ^-m A 2 - ^k 2 ± , (5) 

where we have put p 2 = m 2 , the nucleon mass squared, neglecting the binding. 

The limits of integration in (3) are determined by the condition M 2 > 4m;?, which leads 

to 

x(si + Am 2 ) - -j^—m* A 2 - k\ - Am 2 c > 0. (6) 

J\. X -rl X 

Since k\ > 0, one gets 

x(si + Am 2 ) - -j^—m* A 2 - Am 2 > 0, (7) 

-rl X 

from which one finds the limits of integration in x for the transition A — ► A*: 

x?^ A * <x< x A ^ A * (8) 

where 

x^ A " = —(As- m* A 2 + 4m 2 c ± ^[As - (m^ + 2m c ) 2 ][As - (m\ - 2m c ) 2 ]) (9) 

and s = s\ + Am 2 . The limits of integration in at a given x are determined by (6). 

Using (5) we may pass from the integration variable k\ to \k 2 \. Summing over all states 
of the recoiling nucleus A* we get 



(A) A 2 

a A = / xdx I d\k 2 \a g (x Sl -\k 2 \,k 2 )p(x,\k 2 \), (10) 



where 



*,M,^e(^P) 2 , (id 

\k 2 \min = -r^—X 2 m 2 (12) 

A — x 



and is determined by (9) with m* A put to its minimal value m* A = m A = Am. 

The threshold energy corresponds to = x^ or As = (Am + 2m c ) 2 . In terms of the 
photon energy E we have s\ = 2mE and the threshold energy is found to be 

It steadily falls with A from the nucleon target threshold. With m c = 1.5 GeV we find (in 
GeV): 

Ef = 7.79, Ef = 5.39, Ef = 4.60, = 3.40, E% = 3.02 (14) 
down to the value 2m c = 3 GeV for infinitely heavy nucleus. 
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2.2 High-energy limit 

To interprete p in Eq. (10) it is instructive to study its high-energy limit, which corresponds 
to taking s\ » m 2 and both quantities much greater than the nucleon mass. Assuming 
that the effective values of the gluon virtuality are limited (and small) one then gets for the 
nucleon target (A = 1) 

rl rXSl 

o\ = / xdxa g (xsi) / d\k \p(x, \k |). (15) 

J4m§/si JO 

Here we also neglect the off-mass-shellness of the cross-section off the gluon, considering 
\k 2 \ << Am 2 .. The obtained formula is precisely the standard collinear factorization formula 
with the identification 

f M 2 

xg(x,M 2 )= d\k 2 \xp(x, \k 2 \). (16) 
Jo 

Thus the quantity p(x, \ k 2 \) obviously has a meaning of the double distribution of gluons in 
x and \k 2 \. 



3 The gluon distribution p(x, \k 2 \) 

To find the double distribution of gluons in x and \k 2 \ one may be tempted to use (16) and 
simply differentiate xg(x,M 2 ) in M 2 = \k 2 \. However (16) is only true for x << 1. At 
finite x the derivative dg(x, m 2 )/dM 2 is not positive and cannot be interpreted as the double 
gluonic distribution. 

To avoid this problem, we choose a different, somewhat simplified approach. We assume 
a simple factorizable form for the double density p(x, \ k 2 \) and choose the \k 2 \ dependence in 
accordance with the perturbation theory, with an infrared cutoff in the infrared region: 

P^\ k2 \) = j0A2- (17) 

Function a(x) can be obtained matching (17) with the observed xg(x,M 2 ) at a particular 
point Mq . Since we are interested in the threshold region, we take Mq = 2m c to finally obtain 

/ ,,2 h _ g(x,4mg) 1 , s 

P[X,] IJ " In(4m2/A 2 + 1) |A; 2 | + A 2 - [ ' 

The recipe (18) amounts to taking in (16) p dependent also on M 2 , with the latter dependence 
factorized. Our calculations show that the results are rather weakly dependent of the infrared 
cutoff chosen in the interval 0.4-^0.7 GeV. 

For the proton at x < 1 the gluon distribution g(x,4m 2 ) can be taken from numerous 
existing fits to the experimental data. In our calculations we have used GRV95 LO [4]. For 
the nucleus in the non-cumulative region 

< x < 4 1} (19) 

we use the simplest assumption su(x, Q 2 ) = Agi(x, Q 2 ) neglecting the EMC effect in the first 
approximation. Obviously this not a very satisfactory approximation at x quite close to the 
threshold for the proton target. However our estimates are in any case not justified in this 
region, since, as mentioned in the introduction, multiple gluon exchanges and bound states 
formation begin to play a leading role in the immediate vicinity of the threshold. 

For the nuclei in the cumulative region (outside region (19)) the gluon distribution may 
be estimated using, first, the existing data for the nuclear structure functions in this region 
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and, second, the popular hypothesis that at sufficiently low Q 2 = Qq the sea and gluon distri- 
butions vanish and hadrons become constructed exclusively of valence quarks. Then one can 
find the gluon distribution at a given Q 2 from the standard DGLAP evolution equation, with 
the quark distributions determined from the experimental data on the structure functions at 
x > 1 and evolved back to Q 2 = Qq- In practice we took the initial valence distributions in 
carbon at Q 2 = Qq in the form 

u{x,Ql) = d{x,Ql) = ae~ bx (20) 

and the rest of the distributions equal to zero. Then we calculated the carbon structure 
function at x > 1 and Q 2 in correspondence with the data of [5] and chose the parameters a 
and b to fit the data. With thus chosen a and b we finally calculated the gluon distribution 
in carbon at the scale 4m 2 . Our obtained gluon distributions in carbon for Qq = 0.4 and 0.7 
GeV/c are shown in Fig. 2 for 1 < x < 2. As one observes, the dependence on the choice of 
Qq is very weak in this interval. The slopes result equal to 11.4 (Qq = 0.4 GeV/c) and 11.2 
(Qo = 0.7). 

The distribution for other nuclei was taken from the A-dependence, chosen in accordance 
with the experimental data for hadron production at x > 1 as oc A l+0 - 3x [6]. 

4 Numerical results 

The cross-section (10) involves the photon-gluon fusion cross-section a g off mass shell. The 
integration over the gluon virtuality starts from \k 2 \ ~ m 2 . If one assumes m/M — > then 
the bulk of the contribution will come from the region of small \k 2 \ (with a logarithmic 
precision). In reality m/M is not so small. However, to simplify our calculations, as a first 
approximation, we have taken the photon-gluon cross-section on the mass shell, where it is 
known to be [5] 

2 1 f t2 , \t u AmlM 2 / m 2 M 2 \i 

*W) = I *[; + t + -£-( 1 — srj] < 21 > 

Here e c is the quark charge in units e, u = —M 2 — t and the limits tip are given by 

ii,2 = -y[M± V / M2-4m2] (22) 

We take the strong coupling constant a s = 0.3. 

Our gluon distribution depends on two parameters: the infrared cutoff A in (19) and the 
value of Qq, at which the sea and gluon distributions die out. The order of both is well 
determined, but still one can vary them to some degree. In our calculations we took both A 
and Qq equal to 0.4 or 0.7 GeV/c. 

With these values for the parameters we obtain the cross-sections for charm photoproduc- 
tion on Pb shown in Fig. 3. As one observes, the dependence on both A and Qq is relatively 
weak: in the whole range of their variation the cross-sections change by less than 30%. Fig. 
4 illustrates the A-dependence of the cross-sections (with A = Qq = 0.4 GeV/c). To have 
the idea of the number of nucleons which have to interact together to produce charm at fixed 
energy below threshold we show the limits of intergation xi and X2 in Fig. 5. 

As expected, the cross-sections rapidly fall for energies below threshold. Their energy 
dependence cannot be fit with a simple exponential (in fact they fall faster than the expo- 
nential). As to the aboslute values, for Pb the cross-section fall from ~ 10 pb immediately 
below the threshold down to ~ 1 pb at E = 6 GeV. The A dependence is close to linear. 
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5 Discussion 

We have estimated charm photoproduction rates for nuclear targets below the nucleon target 
threshold. The estimates require knowledge of the gluon distribution in both x and k\ in a 
wide region of the momenta including the confinement region. 

Our estimates were based on a simple factorization assumption and introduction of an 
infrared cutoff. Another approximation has been to take the photon-gluon fusion cross-section 
on its mass-shell. We are of the opinion that this second approximation is not very serious, as 
compared to the first one. In any case it can easily be dropped for the price of considerable 
complication of the calculation. 

Our predictions are infared cutoff dependent. However the cutoff dependence results weak 
for variations of the cutoffs in a reasonable interval. 

In our study we assumed the standard mechanism of charm production via gluon-photon 
fusion (a single gluon exchange between light and heavy quarks). It can be shown that this 
mechanism dominates, provided one is not too close to the threshold (see Appendix). 
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7 Appendix. Multiple gluon exchange 
7.1 Kinematics and phase volume 

To study the relative weight of multiple gluon exchange we consider a simplified picture 
with a mesonic target composed of a light quark and antiquark of mass /i. We neglect the 
binding, so that the meson mass is just 2fi. We shall compare contributions to heavy flavour 
production of the three amplitudes Fig. 6a — a Amplitude a corresponds to a single gluon 
exchange between light and heavy quarks, amplitudes b and c to double gluon exchange. We 
use the light-cone variables and denote ki + = Zip + , pi + = Xip + , i = 1,2. 
The phase volume for the reaction is given by 

16(27r) 8 ziz 2 xix 2 

where d?k\ = dz\d 2 k\_\_ etc and the o-function refers to conservation of the light-cone energy 
(the "-" component of the momentum). Its argument contains the external energy R e = 
2pq+4fi 2 = 2fj,E+4/j, 2 , and the energy of the produced particles R = J2l=i( m c i_i_/ z «+^i_i_/ x i) 
The minimal value of R determining the production threshold occurs at 

Zi = z = c — , Xi = x = ■ — , k i± =p i± = 0, i = 1, 2, (24) 

m c + fx m c + \x 

and is equal to min R = Rq = 2(m c + fi) 2 . 

We shall study our amplitudes near the threshold, so that %i will be small and Zi will 
be close to unity. We put Zi = zq + Q, X{ = xq + i = 1,2 and develop R near the 
threshold keeping terms of the second order in £'s, £'s and transverse momenta. We present 
the difference R e — R in the form R e — Rq = ARq, where dimensionless A measures the 
distance from the threshold and is supposed to be small. Finally we rescale our variables as 
follows 




X^ Rq — 

-, k i± = k i± \/Az R , p i± = p i± \/Ax Ro. (25) 



m 2 
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,w if i \ . m -f X ± 11X2 -r x 2 ) 

V(p 1 ,P2\p' 1 ,p' 2 ) = 4ttq s ^ 2> . (27) 



Obviously new variables with tildes are dimensionless and of the order unity. Using this, one 
finds that approximately £i = — (2 and k 2 = — k\. So the phase volume acquires the form 

dV = y oA 7/2 dz^di^k^dp^dp^ , _ ~ _ ~ 2 _ ~ _ ^ _ x 
(x + eiV2A^)(^o + 6v / 2A^) V ^ 
where V = m 3 M \fz^j2 / (2ir) 8 . 

7.2 Amplitudes 

We use the Coulomb gauge for the interaction between quarks and neglect the contribution 
from the transverse momenta in it. Then the interaction depends only on the scaling variables, 
and for the transition between, say, light quarks p\ + P2 — > p[ + p 2 is given by 

(xi + x' 1 )(x2 + x' 2 ) 
(xi -x[) 

We shall assume that the initial light quarks have their momenta equal to p, so that their 
scaling variable is equal to unity. We omit the the common factor due to their binding into 
the initial target meson. Finally we consider photoproduction, so that q 2 = and choose a 
system in which q + = qj_ = p = 0. 

The amplitude corresponding to Fig. 6a is given by 

A (a) = V(.(P,P\ 2 P-P2,P2)V{q ~ h,2p-p 2 \k 2 ,pi) . , 

{^-{2p-p2) 2 ){ml-{q-k 1 Y) ' 1 ' 

The two interactions near the threshold turn into I2ira s and 2ira s (z2 — z\) where we used 
the fact that x\,X2 « 1 and z\,Z2 ~ 1. For the same reason we find the two denominators 
as /x 2 — (2p — P2) 2 ~ 2 / u 2 ,_i_/ x 2 and m 2 — (q — k\) 2 ~ 2pq. In our dimensionless variables we 
obtain 

^' = -2c lv ^A /f'^'f^ . « = (29) 
fi 2 + 2fi(m c + n)Ap^ ± pq 

The amplitude corresponding to Fig. 66 is 

A (b) = v ((q - h,p\q - k l +p-p 1 ,p l )V(q - ki +p - pi,p\k 2 ,p 2 ) ^ 
(m 2 - (q - ki) 2 )(m 2 - (q - ki + p - pi) 2 ) 

Near the threshold the interactions become — 4ira s and 47tq s The new denominator is 

m 2 c - (q- ki + p-pi) 2 ~m 2 c + 2(m c + ^) /\(^/m~ c ki + v/MPi)!- (31) 
At A << 1 we can drop the second term. So we get 

Ail.) = _„ J_, c 2 = (32) 
m£ pq 

Finally we consider the amplitude of Fig. 6c: 

A (c) = YVM -P + Pi,p\ki,Pi)V(k 2 -p + P2,p\k2,P2) (33 , 
{ml - (ki - p + pi) 2 ){m 2 - (k 2 - p + P2) 2 ) 

Near the threshold both interactions become approximately equal to 47ra s and both denom- 
inators to m 2 . So the amplitude becomes 

A® = 2c 2 ^2c 2 -L (34) 

where we have used that near the threshold pq ~ m 2 . So the amplitudes b and c have the 
same order of magnitude. 
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7.3 Cross-sections 

Now we can pass to our main goal: comparison of contributions of the three amplitudes to 
the total cross-section for heavy flavour production. The first thing to note is that near the 
threshold amplitude a does not interfere with b and c, since since A^ is odd in (i and A^ 
do not depend on £1 at all. Second, since A® and A^ are of the same order and structure 
it is sufficient to compare the contributions of A^ and A^ . Finally, due to the fact that 
xq is small, the magnitude of contributions depends on the relation between A and xq. We 
shall study two limiting cases: A << xq (region A) and A >> xq (region B). 

Region A refers to the production immediately above the threshold. In this case we can 
take xi t 2 — xq. Then we find the contribution of amplitude a to the cross-section as 

= V A 7 / 2 \[2 Cl V2A^ 2 } I (a \ (35) 

where 1^ is a certain integral of the order unity. The contribution of the amplitude b to the 
cross-section will be 

*V = V AVz'[c']>lV, (36) 

where 1^ is another integral of the order unity. The ratio of these two cross-sections will 
have the order 

cr (6) /j, 2 (m c + n) 2 xl 



cjW " m^A ' m 2 c A ~ A 



(37) 



Thus immediately above the threshold the contribution from amplitudes b and c dominate. 
However with the growth of A, in the region x\ « A (and << x to still remain in region 
A) the contribution of amplitudes of b and c become suppressed by factor fi/m c . 

In region B we can approximate xip — £l t 2^/2xoA. To avoid logarithmic divergence in 
£i ; 2 we cutoff the integration region from below at values of the order y/xo/A. We also note 
that the integral over p\± appearing in the contribution of amplitude a is well convergent at 
values of p 2 ± ~ xq/A so that we may neglect p 2 ± in the argument of the 5- function and 
separate the integration over pi± as a factor 

d2pl± * (38) 



+ 2/x(m c + n) Ap 2 ± } 2 2/i 3 (m c + /x) A ' 
We find the cross-section from A^ as 

_ y^^^ -^l— J( «), (39) 

where is an integral of the order ln(xo/A). The cross-section from A^ is found to be 

„<>. = V ^ lc2 -l,M»> (40, 

with jw an integral of the order ln 2 (xo/A). 

The ratio of the two cross-section turns out to be of the same order up to a logarithmic 
factor 

° {b) - ^ + ^ In — ~ In — ~ ^ In — (41) 



cr( a ) m\A xq m 2 A fj, A xq 

and so the contribution of amplitude a clearly dominates in region B, where A >> xq. The 
suppression factor for the contribution of the amplitudes b and c with double gluon exchange 
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is found to be m^/(/j, 2 A). With m c /fi ~ 5 it is of the order 25A. Taking into account 
that double gluon exchange involves a coupling constant at the heavy flavour mass scale will 
add a factor ~3 more. So in the end we find a suppression factor of the order 75A, which 
implies that at a distance of 0.3 GeV from the threshold the contribution of the double gluon 
exchange drops by a factor ~3. 

7.4 Hadrons with more quarks 

Generalization to hadrons with more quarks is straightforward, although quite cumbersome 
due to rapid proliferation of diagrams. However the basic findings remain unchanged. Indeed 
all the difference between amplitides with a single and double gluon exchange between light 
and heavy quarks comes from the fact that a soft propagator of the order x/fij_ in the 
amplitude with a single gluon exchange is substituted by a hard propagator of the order 
l/(pq) in the amplitude with two gluon exchanges. This difference persists for any number 
of quarks in the hadron, although the total number of soft propagator grows. So although 
the overall dependence on A will change (in accordance with the quark counting rules), the 
relation between cross-sections with a single and double gluon exchange will remain the same. 

7.5 Bound states 

One may wonder if the production cross-section is dominated by the formation of final bound 
states, via diagrams as Fig. 7a, which looks as quark rearrangement without any gluon 
exchange [2]. However one has to recall that in the bound state of a light and a heavy quark 
(D-meson) the typical configuration requires pi + /li + = fi/mo and fej+/7j+ = m c /m£>, where 
we neglect the binding taking rri£> = m c + ji. The initial light quarks have however pi+ = p + . 
So for their binding into D-mesons, they have to diminish their longitudinal momenta by at 
least two hard gluon exchanges, as shown in Fig. 7b. But the process in Fig. 7b contributes 
actually a part of the cross-section generated by the amplitude studied in the preceding 
subsection, which corresponds to the immediate binding of the open charm into D-mesons. 
Above the threshold of the open charm production its contribution can only be smaller than 
the total rate of open charm production. True, immediately below this threshold, at distances 
of the order of the binding energy, this mechanism is obviously the only one that contributes, 
in agreement with the estimates above for very small A's. However, as we have seen, with the 
growth of A the strength of multiple interactions between light and heavy quarks necessary 
to produce them in a state appropriate for their binding rapidly goes down. With them goes 
down also the correspondiing part of the cross-section due to immediate binding. 
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8 Figure captions 

Fig. 1. The forward scattering amplitude corresponding to reaction (1). Heavy quarks are 
shown by double lines. 

Fig. 2. The cumulative (x > 1) gluon distributions in carbon at Q 2 = 4m 2 .. The upper 
(lower) curve corresponds to Qq = 0.4(0.7) GeV/c. 

Fig. 3. The charm photoproduction subthreshold cross-sections off Pb for different choice 
of parameters A and Qq. Curves from top to bottom correspond to (A, Qo)=(0-7,0.4),(0.4,0.4), 
(0.7,0.7) and (0.4,0.7) GeV/c. 

Fig. 4. The charm photoproduction subthreshold cross-sections for different targets. 
Curves from bottom to top correspond to A = 12, ,64 and 207. 

Fig. 5. The limits of x-integration for different photon energies and nuclear targets. 
Curves from bottom to top correspond to A =12, 64 and 207. 

Fig. 6. Amplitudes for charm photoproduction off a meson with a single (a) and double 
(6, c) gluon exchanges between light and heavy quarks. The latter are shown with double 
lines. Vertical lines correspond to gluonic exchanges. 

Fig. 7. Amplitudes for the DD photoproduction off a meson. Diagram a is equivalent to 
diagram b, which shows how the produced quarks aquire their momenta appropriate for the 
binding. Notations are as in Fig. 6. 
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